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Abstract. Spatially-resolved high-resolution long-slit spectra of Betelgeuse's circumstellar shell are described for 
a spectral window centered on the 7699A resonance line of neutral potassium. The K I emission from resonance 
fluorescent scattering of photospheric photons which is mapped out to 50 arcsec from the star is approximately 
spherically symmetric with a brightness decreasing as j--2-36±o.03^ where r is the radial distance from the star. 
Our measurements together with the earlier theoretical interpretation by Rodgers & Glassgold suggest that the 
mass loss rate is about 2 xlO~^ Mq . The Kl emission is far from homogeneous: intensity inhomogeneities 
are seen down to the seeing limit of about 1 arcsec and the velocity resolution of about 2 km s~^. There is clear 
evidence for a thin shell of 50 arcsec radius. This is identified with the weaker circumstellar absorption component 
known as S2. Estimates are made of the density of K atoms in this shell (~ 6 x 10~^ cm"''). 
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1. Introduction 

Massive stars may evolve to become red supergiants prior 
to death by supernova explosion. Since mass lost by the 
red supergiant may dictate the time and form of the explo- 
sion and also the nucleosynthetic contribution of the mas- 
sive star to the chemical evolution of the Galaxy, studies 
of mass loss by red supergiants are of great interest. Mass 
loss rates are not yet calculable from an ab initio theory. 
Although measurements of the mass loss rates are avail- 
able from a variety of observations of the gas and dust in 
the circumstellar winds, uncertainties remain large. One 
factor that is often a contributor to the uncertainties is 
an assumption that the stellar surface, circumstellar shell 
and the wind are spherically symmetric. Observations are 
beginning to reveal that supergiants do not recognize fully 
this simplifying assumption. 

Direct observational tests of the assumption are most 
easily executed on the nearest red supergiants. In the case 
of Betelgeuse (a Orionis), the nearest of these stars, de- 
partures from spherical symmetry have been extensively 
noted. At ultraviolet, visible, and near-infrared wave- 
lengths, observations have shown the presence of 'bright 
spots' on or near the stellar surface (cf. Gilliland & Dupree 
1996; Buscher et al. 1990; Young et al. 2000). Such spots 
that change with time are generally attributed to hot ris- 
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ing convective cells, as predicted by Schwarzschild (1975) 
and recently shown by Freytag (2002) in 3D hydrodynam- 
ical simulations (the hot spots appear occasionally over 
downdrafts), but an alternative origin has been proposed 
(Gray 2000). Comparison of ultraviolet and radio bright- 
ness temperatures led Lim et al. (1998) to conclude that 
strong temperature inhomogeneities extend to at least 
several stellar radii above the surface. The measured ra- 
dio brightness temperature at 6cm of about 1400K cor- 
responding to a height of about 7 stellar radii is in good 
agreement with the temperature of about 1200K derived 
by Justtanont et al. (1999) from infrared ground state 
fine structure lines of Sin and Fell. Temperatures de- 
rived from ultraviolet spectra are much higher at these 
heights implying that cold and warm gas coexist at the 
same heights above the stellar surface (Lim et al. 1998). 
Recently, Harper, Brown, & Lim (2001) have constructed 
a semiempirical model for the inner shell extending from 
the photosphere to about 40 stellar radii (about 2 arcsec). 
It is not unlikely that the inhomogeneities detected in ul- 
traviolet and radio observations persist to greater heights 
and into the stellar wind. 

Dust in the shell close to the star has been resolved 
by infrared interferometry and direct imaging, and de- 
tected far from the star by the polarization that scattering 
of dust grains introduces. Infrared emission at about 11 
/im from silicate dust grains comes from an extended re- 
gion of several arcsec with an inner shell radius of 1 to 
2 arcsec (cf. Sloan, Grasdalen, & LeVan 1993; Danchi et 
al. 1994; Bester et al. 1996; Rinehart, Hayward & Houck 
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1998; Sudol et al. 1999). Undoubtedly, dust - cooler grains 

- extends to greater distances but has yet to be resolved 
by far-IR imaging. Extended emission in IRAS images 
was noted by Stencel, Pesce, & Hagcn (1988) and Young, 
Phillips, & Knapp (1993a, 1993b). Warm and cold dust 
is observable by the photospheric light that grains scat- 
ter. This scattered light is separable from light scattered 
by the Earth's atmosphere and the telescope by its polar- 
ization: the polarization vector is tangential to the radius 
vector. McMillan and Tapia (1978) following a suggestion 
by Jura and Jacoby (1976) detected the polarized scat- 
tered light from Betelgeuse's shell. Additional measure- 
ments were made by Le Borgne, Mauron, and Leroy (1986) 
who show that the number density of dust grains declines 
roughly as the inverse-square of the radial distance from 
the star for distances of 10 to 90 arcsec. 

Circumstellar gas as detected by blue-shifted absorp- 
tion lines was observed and analyzed by Weymann (1962) 
who provided indirect evidence that the shell extended to 
great distances from the star. Direct evidence in the case 
of the visual binary Herculis for a very extended shell 
off this M supergiant was presented by Deutsch (1956) 
who found circumstellar lines in the spectrum of the vi- 
sual companion to the supergiant. Weymann's inferences 
about the shell size were compatible with the direct lower 
limit to the outer shell radius from Deutsch's observation. 
Observations of Hi 21 cm emission from Betelgeuse have 
shown that gas extends to about 1 minute of arc from the 
star (Bowers & Knapp 1987). 

Betelgeuse's circumstellar shell was resolved spatially 
by Bernat & Lambert (1975) who detected fluorescent 
emission in the K I 7699A resonance line at angular dis- 
tances of 2 - 4 arcsec from the star. Subsequent obser- 
vations extended the detections of the 7699A emission to 
about 60 arcsec from the star and provided evidence for 
some departures from spherical symmetry and for struc- 
tures in the shell (Bernat & Lambert 1976; Bernat et al. 
1978; Honeycutt et al. 1980; Mauron et al. 1984; Mauron 
1990). Emission in the Na D lines was detected by Mauron 
& Querci (1990), Mauron (1990), and Mauron & Gmlain 
(1995). The most remote (from the star) manifestation of 
the circumstellar shell is the parsec-sized bow shock re- 
ported by Noriega-Crespo et al. (1997) from IRAS images 
at 60 and 100 /im. Emission detected from an arc with a 
mean radius of 6 minutes of arc is attributed to material 
confined by the interstellar medium's ram pressure. 

In this paper, we present and analyse a new series of 
observations of K i 7699A emission obtained with a long 
slit at a resolving power of about 110,000 or 2.6 km s~^ 
with the slit placed at various impact parameters from the 
star and position angles around the star, and at several 
epochs. Our combination of spectral and spatial resolution 
was not achieved by earlier studies of the K i emission. We 
derive afresh the radial dependence of the Ki emission 
and rediscuss the principal conclusion - the mass loss rate 

- that has been previously drawn from the radial depen- 
dence. A novel result of our observations is the discovery 
of clumps in the emission and, hence, in the neutral gas in 



the shell and wind. Some of these inhomogeneities are un- 
resolved both spatially in our seeing-limited observations 
and also kinematically at the limit of about 2 km s^^. 

2. Observations and Reductions 

Betelgeuse was observed in 1994 February, March, and 
November with the W.J. McDonald Observatory's Harlan 
J. Smith 2.7m reflector and its coude spectrograph. The 
latter was used with an echelle grating and an interference 
fllter to separate the order providing the Ki 7699A line. 
The slit length was varied with a majority of the obser- 
vations taken with a length of 10 to 100 arcsec. The slit 
width was almost always set at 0.55 arcsec. Spectra were 
recorded by a Tektronix 512 X 512 CCD over a bandpass 
of about I7A. Standard reduction procedures were applied 
to correct for the bias and pixel-to-pixel variations. The 
wavelength scale was set using Th-Ar lamp spectra. 

During each observing session, several on-star spectra 
were recorded with a very short slit length. Off-star spec- 
tra were then collected with longer slit lengths. In general, 
the star was not included in the slit; inclusion limited the 
exposure time beyond which the starlight saturated a por- 
tion of the CCD chip. Since the fleld at the coude focus ro- 
tates, exposure times were kept short to minimize smear- 
ing of emission inhomogeneities; typically, exposure times 
were 10 minutes or less. For capturing emission far from 
the star, a few exposures as long as 30 minutes were made. 
Relative spectrophotometry was attempted; observations 
were terminated in partly cloudy skies. The stellar ffux 
was measured by opening the slit to its maximum width 
when all the visible image was contained within the slit. 
A guider maintained the star's position relative to the slit 
to an accuracy of about 0.5 arcsec. 

Our reduction techniques were adapted directly from 
those described earlier (Plez & Lambert 1994). In short, 
the spectrum at a given location off the star is presumed 
to be a composite of emission from the gas down the line 
of sight through the shell at that point and starlight scat- 
ter in the Earth's atmosphere and the telescope. We as- 
sume that the latter contribution may be represented by 
a scaled version of an on-star spectrum. Fig. |^ provides an 
example of a raw off-star spectrum and the spectral image 
after subtraction of the photospheric scattered light. For 
positions of the slit near the star (a few arcsec) , the photo- 
spheric scattered contribution largely dominates over the 
light scattered by the shell. This makes the subtraction 
difficult, as the imperfectly removed background spectrum 
shows in Fig. |^. The terrestrial O2 lines appear weakly in 
emission in the image after subtraction. This is due to the 
difference in illumination across the slit when the star is on 
one side of the slit. For spectra taken very near the star, 
we could improve subtraction of the photospheric com- 
ponent by shifting the spectra by less than half a pixel 
in the dispersion direction. Slight changes in airmass be- 
tween the on-star and off-star exposures may also affect 
the subtraction of the O2 lines. 
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Fig. 1. Leftpanel: raw image. Dispersion along the hori- 
zontal axis (l.Skm.s^^/pixcl). Blue to the left. The slit, 
positioned 3" away from the star, spans 100" on the verti- 
cal axis, with a scale of 0.27" /pixel. Right panel: Velocity- 
position map which is the raw image after subtraction of 
the stellar light scattered by the atmosphere and the in- 
strument. Pixels are numbered along both axes. 



The spectrum after subtraction of the scattered light 
is a velocity-position map of the emission. Rows of pixels 
parallel to the dispersion give the spectrum at different 
points in the shell. The radial distance of a point from 
Betelgeuse may be calculated from the impact parameter 
(3 arcsec in this case) and the distance in seconds of arc 
from the point on the slit of closest approach to the star. 
The orientation of the slit was derived from the hour an- 
gle of the observation. A thin uniform shell of outwardly 
flowing and fluorescing K i atoms appears as an ellipse in 
such a velocity-position map. Such an elliptical looking 
feature is evident in Fig. |^. The extent of the ellipse along 
the slit is set by the radius of the shell and the impact 
parameter. The maximum width of the ellipse is directly 
related to the expansion velocity and impact parameter; 
the width vanishes when the impact parameter equals the 
shell's radius. 

Several assumptions underlie our reductions. The 
recorded off-star spectrum results from a convolution of 
the light from points around a position on the slit with the 
seeing profile. A deconvolution of the observed spectra to 
obtain the true spectrum at a point is impossible given the 
varying seeing conditions. We assume that the dominant 
contaminant is scattered starlight with a spectrum rep- 
resented by the on-star spectrum. The on-star spectrum 
is itself, however, not devoid of signatures of circumstel- 
lar K I absorption and emission at 7699 A. In particular, 
there is strong blue-shifted absorption introduced by the 



potassium atoms along the line of sight to the star. At 
some level, this absorption is weakened by emission from 
the immediate environs of the star. Certainly, the angular 
diameter of the stellar photosphere of about 0.05 arcsec 
is considerably smaller than the smallest aperture (0.55 
X 0.55 arcsec) used for the on-star spectra. There is no 
obvious way to retrieve the true stellar spectrum from the 
data. A predicted photospheric profile could be used but 
this would not account for the fact that much of the shell 
absorption probably occurs close to the star. In view of 
the low inferred intensity (relative to the stellar scattered 
light) and the steep gradient of the emitted light from the 
shell, we suggest that our simple reduction technique is 
adequate. 

Our assumption that the circumstellar shell's contribu- 
tion is solely emission in the K i line ignores a contribution 
from the dust grains that scatter photospheric light. This 
scattered light from different locations in the shell will 
be Doppler-shifted by different amounts off the (assumed) 
radially-flowing dust grains, i.e., the scattered light will be 
a smeared version of the photospheric spectrum. Optical 
depth due to dust is small and multiple scattering of pho- 
tons may be neglected. Circumstellar gas is flowing out 
at a velocity of about 10 km s^^. If the dust and gas are 
well coupled, the Doppler shifts imposed on light scat- 
tered by dust grains will amount to a line broadening of 
less than about ±10 km s^^. Since the intrinsic widths 
of the photospheric lines exceed this value, it will be dif- 
ficult to distinguish the dust-scattered contribution from 
that scattered by our atmosphere and telescope. Such a 
distinction is especially hard to achieve with our spectra 
in which the Ki line is the sole strong line. Perhaps, the 
isolation of the dust-scattered component can be made 
using a region richer in deep absorption lines. It may be 
necessary to exploit the polarization of the scattered light 
to effect a clear isolation of the dust-scattered light. We 
assume that the dust makes a negligible contribution to 
the total scattered light. 

Our interpretations are based on the assumption 
that the K i emission results from fluorescence and not 
from radiative recombination or coUisional excitation. 
Calculations using densities and temperatures expected 
of the stellar wind fully support this assumption. 

3. An Asymmetric and Clumpy Shell 

Inspection of the reduced spectra shows that the circum- 
stellar envelope as traced by the Ki emission about the 
star is not completely spherically symmetric. A particu- 
larly striking feature of the velocity-position maps is the 
appearance of clumps at scales extending down to the 
spatial and spectral resolution of our observations. This 
degree of dumpiness has not been revealed by previous 
observations. Sample velocity-position maps are shown in 
Fig. H, H and ^ for slit positions at increasingly larger 
distances from the star illustrate these departures from 
spherical symmetry which we discuss in detail below. 
These Figures also show that the maximum velocity ex- 
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Fig. 2. Velocity-position map. Slit positioned 4" away 
from the star. Dispersion runs horizontally and is 
1.3km.s~"'^/pixel. The slit spans 56" with a pixel size of 
0.27" . The star position is on line 253. Note the very rich 
and clumpy structure, revealing successively larger incom- 
plete emission shells. At this position on the sky, around 
column 260, there is more red-shifted than blue-shifted 
emission. Pixels are numbered along both axes. 



tent of a velocity-position map is not greatly different for 
small and large impact parameters. 

3.1. Radial Dependence of the Ki Emission and Mass 
Loss Rate 

By azimuthal averaging to the maximum extent possible 
for a given sequence of observations, it is possible to deter- 
mine the mean radial dependence of the emission intensity 
integrated over the line profile. This was done only for ex- 
posures obtained under good photometric conditions. The 
intensity of the scattered light at the radial distance r is 
written as Iscatt{r). We scale the intensity to the mea- 
sured stellar flux /, measured with the very wide slit. /, 
was measured by integrating the on-star spectrum flux 
over the whole spectral interval (512 pixels corresponding 




Fig. 3. Two observations at the same orientation and 
slit position (13" from Betelgeuse). Velocity scale as 
in Fig. 0: 0.27"/pixel along the slit(vertical axis), and 
1. 3km. s /pixel on the dispersion axis (horizontal). The 
star position corresponds to line 190. Left panel: March 
22, 1994 with 1.2" seeing. Right panel: Nov 28, 1994, with 
seeing around 3" . Note the great similarity of the images. 
Most differences are probably due to seeing degradation. 
Pixels are numbered along both axes. 



to 17.24A). This total flux was then divided by the wave- 
length range and the integration time. The resulting 
is thus expressed in counts per A per second. Note that 
Mauron et al. (1984) used a similar definition in averag- 
ing the flux over 20A, while Mauron (1990), Honeycutt 
et al. (1980), and Bernat et al. (1978) used the photo- 
spheric flux in lA centered on 7699A. There are a number 
of sources of uncertainty affecting the determination of 
the photospheric flux and in the Iscattir)/ 1* ratio, mostly 
because of possible loss of photons intercepted by the slit 
jaws, and changes of transparency and seeing conditions 
between exposures. The derived normalized intensity of 
the K I emission, Iscatt (r) //* , is shown in Fig. || for im- 
pact parameters R out to 50 arcsec from the star. The 
mean profile is consistent with a power law of slope -2.36 
± 0.03 (least-square fit to the data). Perhaps, one should 
fit the data with a two component model: a slope of -2.59 
± 0.2 out to 7 arcsec and of 2.55 ± 0.07 beyond 10 arcsec 
but brighter by a factor of 1.4. 

Our measured slope is equal within the uncertain- 
ties to that reported earlier by Honeycutt et al. (1980) 
who gave the slope as -2.65 ± 0.2 from average measure- 
ments in the four principal quadrants at radial distances 
of about 10 to 50 arcsec. Estimates of the intensity of 
the scattered light are similarly in good agreement: we 
obtain Iscatt/ 1* — 10~^ ° at 20 arcsec from the star, and 
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Fig. 4. The S2 shell. Sht at 33" from the star. Scale as 
in Fig. H Star position corresponds to line 261. Pixels are 
numbered along both axes. 

Honeycutt et al. reported a value of 10^^'^. The agreement 
is less satisfactory with Mauron et al. (1984) who reported 
a slope of -3.5 ±0.8 from measurements spanning distances 
of 20 to 55 arcsec but the intensity at 20 arcsec (10~^'*) is 
consistent with other measurements. Mauron's (1990 - see 
also Mauron & Querci 1990) measurements at 5 to 8 arc- 
sec are similarly a little higher in intensity than our mean 
values. Although the differences could be due to changes 
in the efficiency of fluorescence, the density of potassium 
atoms in the shell, and other physical factors, we think it 
more likely that they reflect observational errors encurred 
in the surface photometry. As noted by Mauron (1990), 
earlier results by Bernat et al. (1978) gave much higher 
intensities for the emission close to the star than are now 
reported: e.g., Bernat et al. found Iscatt/h — 10"^ '^ at 




0.5 1 1.5 

log [(9{arcsec)] 



Fig. 5. Average scattered intensity profile. Iscatt/I* is the 
intensity scattered by the shell, normalized to the photo- 
spheric flux (see text for details). Eleven profiles from ex- 
posures at various distances from the star were averaged. 
In each individual image, the intensity was integrated in 
1" bins along the slit. Error bars are one standard devi- 
ation. The best linear fit to the data out to 50 arcsec is 
superimposed (slope=-2.36). The position of the 50 arcsec 
shell is indicated by the tick mark. Note the strong drop 
in intensity beyond that distance to the star. 

about 3 arcsec but our value is lO^'^'^. This large differ- 
ence is unlikely to be due to changing conditions at the 
stellar surface or in the shell but reflects the difficulty of 
making photometric measurements of a faint light source 
close to a bright star with a detector (a vidicon) with a 
non-linear response. 

A calculation of the stellar mass loss rate from the 
radial dependence of the 7699A emission has been dis- 
cussed previously, notably by Rodgers & Glassgold (1991). 
Although the observations may provide reliable estimates 
of the run of the density of potassium atoms, an estimate 
of the mass loss rate calls for knowledge the degree of ion- 
ization of potassium, the depletion of potassium onto dust 
grains, excitation of the atoms (fluorescence vs radiative 
recombination, for example) in the shell, and radiative 
transfer through the shell. We assume the shell to be op- 
tically thin to the resonance photons, a conclusion that is 
not contradicted by the ratio of emission in the two reso- 
nance lines (Mauron & Querci 1990). Here, we adopt the 
recipes provided by Rodgers & Glassgold and adjust their 
derived mass loss rate to our measured intensities. 

The slope of the Iscauif)/!* versus r relation is, if 
spherical symmetry is assumed, directly related to the ra- 
dial dependence of the number density of neutral potas- 
sium atoms. Additional assumptions include uniform illu- 
mination of the atoms by the photosphere, and a shell that 
is optically thin to the 7699A photons. It is then readily 
shown that a relation N{r) oc provides an intensity 
dependence Iscatt/I* oc r^(^+^) (Bernat 1976; Gustafsson 
et al. 1997). If neutral potassium atoms are conserved in 
a radial outflow, their density dependence follows 9 — 2 
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and I scan /I* follows a power law of slope 3. Conversely, 
the observed slope of -2.36 indicates a density dependence 
with a slope of -1.36. Our observed slope is in good agree- 
ment with the value of -2.33 corresponding to the model 
by Rodgcrs & Glassgold. The difference between their pre- 
diction and the slope of -3 predicted for a uniform shell 
with conservation of potassium atoms is due to their in- 
troduction of ionization of potassium atoms in the shell by 
chromospheric ultraviolet photons. Glassgold & Huggins 
(1986) and Rodgers & Glassgold (1991) discuss in detail 
the run of the neutral potassium fraction. The ionizing 
chromospheric radiation field is not much shielded in the 
envelope and decreases as r^^. The electron fraction, set 
constant in Glassgold & Huggins model, increases slightly 
outwards in Rodgers & Glassgold's more detailed model 
resulting in a gradient of the electron density which is shal- 
lower than r^^ . The recombination coefficient increases 
with distance as a result of the temperature decrease. The 
neutral fraction thus increases with distance from the star 
in the inner shell. At large distances out from the star 
where ionization of potassium atoms is controlled by the 
interstellar radiation field, the slope of Iscattif)/ 1* must 
steepen. This steepening is visible beyond about 30 arcsec 
in our data on Fig. ||, although at first sight masked by 
the detached shell S2 (see Sec. 3.2.2). The slope increase 
(to about -4.2) is consistent with Rogers and Glassgold 
(1991) model. 

Following Gustafsson et al. (1997) we write: 
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where /3 is the angular distance from the star in arc- 
sec, / is the oscillator strength of the line, Ngcatt is the 
number density of scattering atoms in cm~'^, d is the dis- 
tance to the star in parsec, and z is the distance along 
the line of sight in cm. Assuming a power law of the 
form: Nscatt = Nsisi^^)^ , "^ith TV^is the number den- 
sity at r = lO^^cm, the integral is readily calculated and 
yields: iVsi5(10^^)^y(22|f5)(e+i)^ ^ith Y a constant de- 
pending on the value of 9. From our observations we derive 
^^s^fi^ = "-"^V™"" ■ Adopting a distance of 140 pc for 
Betelgeuse from its Hipparcos parallax of 7.63 ± 1.64 mas, 
we find the number density of scattering atoms to be 
iV(Kl) = 6.0 xlO-4 X (IQiVO^-^^ cm-3.[] 

The column density of neutral K atoms along the line 
of sight to the star was derived by Bernat (1977) from the 
7699 A circumstellar absorption line: 7V(Kl) ~ 4 x 10^^ 
cm^^. Integrating our relation for the radial dependence, 
Bernat 's column density is reached by integrating from 



^ Gray (2000) has remarked that the accuracy of the 
Hipparcos parallax may be seriously underestimated because it 
assumes that the centroid of the stellar image remained fixed 
relative to the stellar disk over the period of observations; a 
bright spot midway between disk center and the stellar limb 
and contributing 25% of the light could cause the centroid to 
shift by half the parallax. 



infinity to an inner radius Rmner = 9 x IO^^cto ~ 2Rstar- 
This is a plausible value for the inner radius. To halve 
the column density requires Rinner — 13i?star- To double 
for Rinner < Rstar, which makcs no sense. In short, our 
derived radial dependence for the density of potassium 
atoms is consistent with the column density derived from 
the absorption line. 

For this reason and because our observed slope for 
Iscatt{r)/ 1* is equal to that predicted by Rodgers & 
Glassgold, we use the ratio of our observed to their pre- 
dicted number density versus radial distance directly to 
correct their mass loss rate. This corresponds to a reduc- 
tion of their mass loss rate from Af = 4 x 10~® Mq 
to 2.4 x 10~^ Mq y~^. The revision is dependent on our 
value of Iscatt 1 1* which rests on an average of the stellar 
photospheric flux averaged over 17 A around the potas- 
sium line. The stellar flux intercepted by the outflowing 
gas of the shell is that on the blue wing of the Ki line. 
As already pointed out, we do not have access to the true 
photospheric spectrum, and the K i line profile can only 
be guessed at. A rather conservative estimate of a factor of 
2 lower flux at line center leads to a mass-loss rate larger 
by a factor of two. A distance of 200 pc was assumed by 
Rodgers & Glassgold and Glassgold & Huggins. To correct 
the mass loss rate to the Hipparcos distance oi d = 140 
pc, we use the scaling given by Mauron & Caux (1992), 
i.e., the Ki brightness at a given angular distance from the 
star is approximately proportional to /d. Adoption of 
d — 140 pc instead of 200 pc implies a reduction of M 
by 85 % or Af = 2 x IQ-^Mq y-^. Glassgold & Huggins 
-2 ^arn that the accuracy of an M derived from K i emission 
is low ( "a factor of some 3 - 5" ) . Although the subsequent 
analysis of the temperature profile (Rodgers & Glassgold 
1991) may have reduced the uncertainty, a factor of two 
uncertainty seems likely to be optimistic. The derived rate 
is dependent too on the assumption of a spherically sym- 
metric shell, an assumption probed observationally in the 
next section. 

This uncertainty could in principle be reduced by ob- 
serving other emission lines. In light of the higher abun- 
dance of Na and using the atomic concentrations predicted 
by Rodgers & Glassgold's (1991), the Na D emission is ex- 
pected to be stronger than that in the K i 7699 and 7665 A 
lines. Mauron & Guilain (1995) find that the Na D emis- 
sion is considerably fainter than expected although other 
mass losing M stars do show Na D emission stronger than 
the Ki emission. These authors suggest that interstellar 
absorption is "the simplest explanation, though not com- 
pletely convincing" . Fortunately, interstellar absorption in 
the K I lines is observed to be considerably weaker than in 
the NaD lines. 

Betelgeuse has been subjected to numerous investiga- 
tions of its mass loss rate using a variety of spectroscopic 
indicators. The Ki-based rate may be compared with a 
couple of determinations using tracers of a major con- 
stituent of the shell. Bowers & Knapp (1987) detected 21 
cm Hi emission to derive M ~ 1.1 x 10~® Mq if d 
= 140 pc. Huggins et al. (1994) detected the Ci 609 /^m 
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fine structure line and inferred a mass loss rate also of 
M ^ 1 X 10"^ Mq y-i for d = 140 pc. The close agree- 
ment between these results and ours is encouraging but 
must be largely fortuitous. More recently, Harper et al. 
(2001) estimated M = 3.1 ± 1.3 x lO-^M© y-^ from ab- 
sorption features in ultraviolet Fe ii lines. 

3.2. Departures from Spherical Symmetry 

Azimuthal averaging produces a smooth run of the emis- 
sion line's flux with radial distance but masks a fascinat- 
ing structure of incomplete shells and clumps. This struc- 
ture is illustrated in Fig. |[ and ^ and additionally in 
Fig. H, and 0. Several features of the structure are seen 
by inspection: the appearance of thin incomplete shells, 
intense clumps of small spatial extent and velocity disper- 
sion, and lines of sight - particularly far from the star and 
inside the outermost detected shell - from which no 7699A 
emission is detected, and a seemingly sharp boundary to 
the outermost shell. 

The following relations are helpful in discussing the 
properties of the circumstellar shell. At the assumed dis- 
tance of 140 pc, 1 second of arc corresponds to a projected 
distance of 140 AU or a light travel time of about 19 hours. 
The stellar angular diameter is taken to be 0.044 arcsec 
(Dyck et al. 1998) or a stellar radius of 3.1 AU. At an 
expansion velocity of 15 km s^^, gas travels 3.2 AU in 1 
year or an angular distance of just 0.02 arcsec. 

3.2.1. Temporal Dependence? 

Our observations were obtained over several months in 
1994. In different observing runs, several successful at- 
tempts were made to acquire a series of spectra with the 
slit set relative to the star at the same positions and ori- 
entations. A goal of these observations was to search for 
changes in the velocity-position maps. 

Expected time scales for changes to occur in the wind 
itself are much longer than a few months (see below) . More 
rapid changes may result from the rise and fall of bright (or 
dark) spots on the stellar surface. High angular resolution 
images of Betelgeuse's disk show a few (1 - 3) bright spots 
contributing 10 - 20 per cent of the light and changing on 
a timescale of weeks to months (Wilson, Dhillon, & Haniff 
1997). With a light travel time of about a month to a ra- 
dial distance of about 30" , intensity fluctuations of (say) 
20 % evolving on a time scale of a month are not going 
to produce dramatic changes in the fluorescent emission. 
Were a spectacularly bright spot to occur, a large inten- 
sity contrast would be seen in the hemisphere facing the 
spot. The appearance of shells and especially the presence 
of clumps cannot be traced to spots on the illuminating 
photosphere. 

Changes in the density of potassium atoms depend 
on the rates of ionization and recombination. With rates 
taken from Glassgold & Huggins (1986), we find the 
timescale for photoionization by the stellar radiation field 
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Fig. 6. Two observations at the same orientation and slit 
position (21" away from the star). The star corresponds 
to line 191. Exposure time 30 min. Left panel: March 23, 
1994 with 1.7" seeing. Right panel: November 26, 1994 
with seeing around 1.2". The general appearance of the 2 
frames is similar, but more differences are visible than on 
Fig. ^. The central condensation around pixel (5,210) is 
more conspicuous in the left-hand frame, and differences 
in the clump structure are obvious in the region around 
(5,270). 

to be tphoto ~ 0.16*^ yr where 9 is the distance from the 
star expressed in arcsec. The recombination timescale is 
necessarily much longer: tree ^ 306*^ yr (or ^^^^ with the 
temperature dependence of the recombination coefficient 
taken into account). Except perhaps very close to the star, 
these estimates show that changes in K i emission occur- 
ring on a time scale of a few months are unlikely to be 
traceable to variations in the photoionizations by the stel- 
lar radiation field. 

Fig. ^ shows pairs of velocity-position maps taken at 
the same slit position relative to the star but at 8 months 
apart. By inspection, there is considerable similarity be- 
tween the March and November observations for which the 
impact parameter of the slit was 13 arcsec. Gross struc- 
tures may be deemed identical. Smaller structures (here 
dubbed 'clumps') show no discernible change in location. 
Differences in shape of the smaller structures may be at- 
tributed to differences in the seeing: the seeing was about 
1.2 arcsec for the March observation but about 3 arcsec 
in November. 

Fig. ^ and show similar pairs of velocity-position 
maps from the March and November runs but with an im- 
pact parameter of 21 arcsec. The exposure in Fig. was 
obtained immediately following that in Fig. ^. In these 
figures, there is evidence for changes in the smaller struc- 
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Fig. 7. Same as Fig. ^. These spectra were acquired di- 
rectly following the ones of Fig. at the same slit posi- 
tion. While there is considerable similarity between this 
and Fig. ^, there are some differences. These illustrate the 
changes in appearance revealed by the rotation of the slit 
on the sky. Comparison of the right-hand panels of the 
two figures shows differences in the regions around pixels 
(55,80), (55,250), and (50, 330). 



tures. Indeed, the March pair of maps clearly differ in the 
inner parts from the same area in the November maps. 
Close inspection shows a few differences for clumps in suc- 
cessive maps in March and in November. This is not sur- 
prising given that successive maps owing to image rotation 
at the coude focus do not sample exactly the same por- 
tion of the shell. The telescope was guided on the star and 
the slit rotates about the star. Changes in position on the 
sky recorded by a particular pixel can exceed the angular 
size of some clumps, and, hence, these may appear clearly 
in one image but not in the next. This instrumental effect 
increases with increasing impact parameter which may ex- 
plain in part why changes are more evident in Fig. ^ and 
I than in Fig. | 

3.2.2. A Distant Spherical Shell 

A striking feature of the velocity-position maps with im- 
pact parameters greater than about 20 arcsec is the ap- 
pearance of an elliptical shape for the outermost emission 
(Fig. 1^. Emission is very weak, possibly absent, beyond 
about 1 minute of arc from the star. EUipticity of isophotes 
in the velocity-position map corresponds to a spherically 
symmetric expanding shell. Geometry of the slit position 
(at 33 arcsec from the star) and the semi-major axis of 
the ellipse show that the shell radius is 50 ± 2 arcsec. The 
maximum velocity width, the ellipse's semiminor axis, cor- 



responds to an expansion (radial) velocity of 18 ± 2km 
s~^. Note that the shell is clearly seen in Fig. || as the 
plateau followed by an intensity drop at about 50 arcsec. 
At that distance from the star, the mean Ki density is 
rapidly decreasing with distance, due to the interstellar 
UV field (see Fig. |5|and Sec. 3.1), and the shell represents 
a considerable density enhancement. 

The thickness of the shell is 1 to 2 arcsec in position on 
the sky and is unresolved in velocity, i.e. less than 2.6 km 
s~^. Non-uniform emission over the spherical shell could 
result from one or more of several factors: shell ejection 
occurred preferentially over one hemisphere of the star; 
the excitation of the shell is non-uniform owing to a bright 
spot on the stellar disk; the ionisation of potassium atoms 
is greater over one hemisphere due to a chromospheric 
enhancement; the shell is plowing into the irregular local 
interstellar medium and deceleration enhances the local 
density non-uniformly. 

Circunistellar absorption imposed upon the photo- 
spheric spectrum is composed of two components in 
atomic (e.g., Ki 7699 A as shown by Bernat & Lambert 
[1975] and Goldberg et al. [1975]) and the CO lines by 
Bernat et al. (1979) who labelled the components SI and 
S2 where the former has a higher excitation temperature 
(200K vs 70K), substantially larger CO and Ki column 
densities and lower expansion velocities. The expansion 
velocity deduced from the S2 absorption lines (20.2 km 
s~^ from Goldberg et al., and 13 km s~^ from Bernat et 
al.) are consistent with our measurement. Note that the 
absorption line measurements are made relative to pho- 
tospheric lines but the photospheric velocity is variable. 
Our measurement is a true expansion velocity, relative to 
the center-of-mass of the star, independent of the stellar 
velocity, and of photospheric velocity variations. 

Another notable difference between the SI and S2 com- 
ponents is the very narrow line width of the latter: a 
Doppler width of 1 km s~^ for S2 but 4 km s^^ for SI 
from visual lines (Bernat 1977). The Doppler width for S2 
corresponds to a full width at half intensity of 1.7 km s~^, 
an estimate consistent with the upper limit obtained from 
our observation of the shell emission. Bernat et al. from 
a comparison of the coUisional and radiative rates for CO 
excitation inferred that the former dominate and, hence, 
their excitation temperature was approximately equiva- 
lent to the gas kinetic temperature. Reference to a dust 
shell model (Tsuji 1979) showed a gas temperature of 70 K 
was achieved at about 55 arcsec from the star, a prediction 
quite consistent with our observed shell radius. 

Judged by velocity of expansion, line width, and shell 
radius as measured or estimated from the CO lines and our 
7699A velocity-position maps, we identify the circumstel- 
lar gas contributing the S2 absorption lines with the out- 
ermost thin shell seen in Fig. |4[ This identification allows 
interesting further deductions to be made. In particular, 
the volume density of CO molecules and neutral K atoms 
may be estimated from their column density and the shell 
thickness. The equivalent width of the K i 7699A S2 line 
estimated from Goldberg et al. (1975, their Figure la) 
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gives a column density N(K i) ~ 1.3 x 10^^ cm^^ on the as- 
sumption that the hne is unsaturated. After correction for 
seeing, the shell's thickness may be about 1 second of arc. 
At d = 140 pc, a shell thickness of 1 second of arc corre- 
sponds to 2 X 10^*^ cm. Then, the K i density in the shell is 
about 6 xlO~^ cm~'^. In contrast, Rodgers & Glassgold's 
(1991) model predicted a density of n(K i)~ 2 x 10~^ cm~'^ 
at the shell's distance. This difference is not surprising as 
the shell by inspection of Fig. ^ is evidently a considerable 
density enhancement above the gas at smaller and greater 
distances In addition, the intensity variation around the 
shell is obvious implying that the cut through the shell 
to the star that gives the K i absorption line may not be 
representative of mean conditions in the shell. 

If the shell at 55 arcsec is identified with Bernat et al.'s 
S2 shell, their measured column density (1.2 xlO^^ cm^) 
translates to a volume density n(CO)~ 6 cm^^. The ratio 
of the CO and K i column densities suggests a density ratio 
n(CO)/n(Ki)~ 9 x 10'*. The elemental abundance ratio 
in the photosphere is probably a little less than the solar 
ratio of 2700 on account of a lower C abundance resulting 
from the first dredge-up. If C has been reduced by about 
50%, and K is not depleted onto grains in the wind, the 
expected elemental ratio in the wind is 1350. The higher 
ratio in the S2 shell implies, as expected, that K is largely 
ionized: the fraction of neutral K atoms is j/k — 0.014 
provided that most of the C is in CO. Huggins et al.'s 
(1994) prediction of ionization of C and association into 
CO through the shell suggests CO accounts for about 20 
% of the carbon. Then, t/k 0.003. Rodgers & Glassgold 
(1991) 'standard' model predicts t/k — 0.01 at the radial 
distance equivalent to 50 arcsec at d= 140 pc. This is fair 
agreement with observation given the various (uncertain) 
ingredients that enter into both estimates. The mass in 
the shell amounts to about 10~^ M0, at a distance of 140 
pc and assuming ?/k — 0.01. 

As noted above the S2 shell is not a consequence of 
a varying photospheric illumination; the light travel time 
across the shell is only about a day, and the shell is only 
about 50 light days from the star. In addition, the shell 
seen in absorption in the 1970s (Goldberg et al. 1975; 
Bernat et al. 1979) is plausibly identified with the Ki 
emission shell seen in 1994. The simplest explanation for 
the shell may be that it is mass ejected at a substantially 
higher than average rate at an earlier time. If the expan- 
sion velocity has been constant, gas travelling at 18 km 
s~* and now at 50 arcsec from the star was ejected 1900 
years ago {d — 140 pc is assumed). A shell thickness of 
1 second of arc implies ejection lasted about 40 yr but 
this is a guess because the velocity of ejection presum- 
ably varied too. The shell as it appears in Fig. ^ is not 
uniformly bright. It is unlikely that the observed contrast 
around the shell is due to the presence of bright spot on 
the disk at the time of the observations. A possible expla- 
nation is that mass ejection occurred preferentially above 
a spot 1900 years ago. Stellar rotation (the rotation pe- 
riod of Betelgeuse is at least a few years) , could have led 
to material being ejected into the shell at a variable rate 



over the lifetime of the spot. There appear to be at least 
two other shells of moderate completeness which suggests 
that significant increases in mass loss may occur every 
several hundred years. Multiple shells were also detected 
around fi Cep by Mauron (1997), with a characteristic 
time scale between two successive increases in mass loss 
of about one thousand years. Detached shells are also ob- 
served around AGB stars, where they are attributed to a 
Helium shell flash that induces both a mass loss increase 
and a two-wind interaction due to the increased outflow 
velocity (Steffen and Schonberner 2000). Clearly, an alter- 
native mechanism must be operating in supergiant stars 
like Betelgeuse. 

Future observations of the K i emission from the shell 
should be made to search for changes in the S2 shell. 
Measuring the S2's expansion will be a challenge: the 
proper motion expected is a mere 0.03 arcsec per year. 

3.2.3. Clumps 

A striking feature of our velocity-position maps, especially 
those acquired in the best seeing is the appearance of 
clumps or knots of emission; Fig. ^, |^ and show clumps 
extending in size down to the velocity and spatial resolu- 
tion of the map. There are locations where clumps appear 
connected suggesting they are condensations in a sheet or 
shell. In general, the collection of clumps in many velocity- 
position maps accounts for more than half of the emit- 
ted flux. The impression of elliptical structures is largely 
gained from gas at large angular distances from the star. 
Many of our exposures, especially those taken with small 
impact parameters, were too short to reveal faint distant 
structures. 

These are the first observations of such clumps in 
Betelgeuse's wind. Mauron (1997) from long-slit spec- 
tra showing K i emission from fi Cephei reports 'inhomo- 
geneities' with an apparent size of 1" to 3" when the an- 
gular resolution was about 0.6". In velocity space, /x Cep's 
clumps were unresolved at the moderate resolution of 40 
km s^^ of these observations. Since n Cep is about 5 times 
more distant than Betelgeuse, its smallest clumps would 
appear to be 5 to 10 larger than those we have detected. 

There is a hint that the prevalence of clumps is higher 
close to the star than at the largest distances probed 
by our spectra. This encourages the speculation that the 
clumps' origin may be traceable through the inhomoge- 
neous inner shell (Lim et al. 1998) and chromosphere 
to the stellar surface. The greater visibility of shells at 
large distances may be due to the dissolution and merger 
of clumps at a common distance. Although clumps are 
seen by fiuorescent emission, their shape and size are 
attributable to an enhanced concentration of potassium 
atoms in the wind and not to intensity variations over the 
surface of the illuminating photosphere; a bright photo- 
spheric spot will enhance fluorescence over approximately 
2tt of solid angle not a mere few arcsec. The clumps are 
probably not regions of reduced ionization of potassium 
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atoms because the characteristic timescale for photoion- 
ization is shorter that the expansion timescale (the time 
required to reach a given angular distance at a constant 
expansion velocity): it is a factor of 100 faster close to the 
star and about a factor of 10 at 50 arcsec from the star. 

Supposing the clumps represent coherent structures 
related to the process of mass loss, a change of clump 
characteristics with radial distance provides information 
about the clumps' internal structure. Consider a clump 
moving radially outwards at velocity Vexp with an internal 
velocity dispersion Vdiap- Its radial size cp increases linearly 
with radial distance out from the star i.e., A(/)/A<l> = 
Vdisp/vexp- Identifying Udisp with the sound speed, we esti- 
mate Vdisp = 1 km s~^. Since Vexp — 15 km s~^, a clump of 
about 2 arcsec at 10 arcsec from the star would expand to 
about 5 arcsec at a radial distance of 50 arcsec. In the in- 
ner shell, the inter-clump distances are a few arcsec. Then, 
if dispersion of the clump perpendicular to the radial flow 
occurs at about 1 km s~^, the clumps at a common dis- 
tance may well merge after a few tens of arcsec. 

Our observations are consistent with this result in that 
the maximum velocity extent of a velocity-position map 
is not greatly different for small and large impact param- 
eters. 

4. Concluding Remarks 

Our study of the Ki emission from Betelgeuse's circum- 
stellar shell reveals new aspects of the shell and its potas- 
sium atoms. Spherical symmetry is no more than a rough 
approximation to the distribution of potassium atoms but 
the radial distribution of the brightness closely follows the 
predictions for a model proposed by Rodgcrs & Glassgold 
(1991) implying a mass loss rate of about 2 xlO~^ Mq 
y~^. But the striking result of our study is the presence 
of clumps of potassium atoms with the smallest chimps 
being unresolved spatially and kinematically. The clumps 
likely originate from cool structures close to the star. 
Observations suggest that they disappear at large dis- 
tances from the star. 

Observations at a site of excellent seeing are to be 
sought in order to trace the clumps back closer to the 
star. Other supergiants and giants are observable in both 
Na D and K i emission, as published reports demonstrate, 
but the occurrence of clumps in the winds from such stars 
is as yet unknown. To explore the innermost regions of 
these shells will require the spatial resolution of a space, or 
a diffraction limited groimd-bascd telescope. Betelgeuse's 
chromosphere has already been explored in the ultravi- 
olet continuum and the Mgll h and k lines (Uitenbroek, 
Duprcc, & Gilliland 1998). Extension of this work to other 
lines, and other stars would be rewarding. 

Acknowledgements. We thank N. Mauron and E. Josselin for 
inspiring discussions of the Betelgeuse shell. 

References 

Bernat, A. P. 1976, Ph.D. thesis. University of Texas at Austin 



Bernat, A.P. 1977, ApJ, 213, 756 

Bernat, A.P., Hall, D.N.B., Hinkle, K.H., & Ridgway, S.T. 

1979, ApJ, 235, L135 
Bernat, A. P., Honeycutt, R.K., Kephart, J.E., Gow, C.E., 

Sandford, II, M.T., & Lambert, D.L. 1978, ApJ, 219, 532 
Bernat, A.P., & Lambert, D.L. 1975, ApJ, 201, L153 
Bernat, A.P., & Lambert, D.L. 1976, ApJ, 210, 395 
Bester, M., Danchi, W.C., Hale, D., Townes, C.H., Dcgiacomi, 

C.G., Mekarnia, D., & Geballe, T.R. 1996, ApJ, 463, 336 
Bowers, P.P., & Knapp. G.R. 1987, ApJ, 315, 305 
Buscher, D.F., Haniff, c.A., Baldwin, J.E., & Warner, P.J. 

1990, MNRAS, 245, 7 
Danchi, W.C., Bester, M., Degiacomi, C.G., Greenhill, L.J., & 

Townes, C.H. 1994, AJ, 107, 1469 
Deutsch, A.J. 1956, ApJ, 123, 310 

Dyck, H.M., van Belle, G.T., & Thomson, R.R. 1998, AJ 116, 
981 

Freytag, B. 2002, in Cool Stars, Stellar Systems, and the Sun, 

12th Cambridge Workshop, in press. 
Gilliland, R.L., & Dupree, A.K. 1996, ApJ, 463, L29 
Glassgold, A.E., & Huggins, P.J. 1986, ApJ 306, 605 
Goldberg, L., Ramsey, L., Testerman, L., & Carbon, D. 1975, 

ApJ, 199, 427 
Gray, D.F. 2000, ApJ, 532, 487 

Gustafsson, B., Eriksson, K., Kiselman, D., Olander, N., & 

Olofsson, H. 1997, A&A, 318, 52 
Harper, G.M., Brown, A., & Lim, J. 2001, ApJ, 551, 1073 
Honcycutt,R.K., Bernat, A. P., Kephart, J.E., Gow, C.E., 

Sandford, II, M.T., & Lambert, D.L. 1980, ApJ, 239, 565 
Huggins, P.J. 1987, ApJ, 313, 400 

Huggins, P.J., Bachiller, R., Cox, P., & Forveille, T. 1994, ApJ, 
424, L127 

Jura, M., & Jacoby, G.H. 1976, Ap Lett, 18, 5 

Justtanont, K., Tielens, A.G.G.M., de Jong, T., Cami, J., 

Waters, L.B.F.M., & Yamamura, I. 1999, A&A, 345, 605 
Le Borgne, J.F., Mauron, N., & Leroy, J.L. 1986, A&A, 168, 

211 

Lim, j., Carilh, C.L., White, S.M., Beasley, A., Marson, R.G. 

1998, Nature, 392, 575 
Mauron, N. 1990, A&A, 227, 141 
Mauron, N. 1997, A&A, 326, 300 
Mauron, N., & Caux, E. 1992, A&A, 265, 711 
Mauron, N., Fort, B., Querci, F., Dreux, M., Fauconnier, T., 

& Lamy, P. 1984, A&A, 130, 341 
Mauron, N., & Guilain, Ch. 1995, A&A, 298, 869 
Mauron, N., & Querci, F. 1990, A&AS, 86, 513 
McMillan, R., & Tapia, S. 1978, ApJ, 226, L87 
Noriega- Crespo, A., van Buren, D., Cao, Y., & Dgani, R. 1997, 

AJ, 114, 837 
Plez, B., & Lambert, D.L. 1994, ApJ, 425, LlOl 
Rinehart, S.A., Hayward, T.L., & Houck, J.R. 1998, ApJ, 503, 

415 

Rodgers, B., & Glassgold, A.E. 1991, ApJ, 382, 606 
Schwarzschild, M. 1975, ApJ, 195, 137 

Sloan, G.C., Grasdalen. G.L., & LeVan, P.D. 1993, ApJ, 404, 
328 

Steffen, M., & Schonberner, D. 2000, A&A, 357, 180 
Stcncel, R.E., Pesce, J.E., & Hagen, W. 1988, AJ, 95, 141 
Sudol, J.J., Dyck, H.M., Stencel, R.E., Klebe, D.L, & Creech- 

Eakman, M.J. 1999, AJ, 117, 1609 
Uitenbroek, H., Dupree, A.K., & Gilliland, R.L. 1998, AJ, 116, 

2501 

Wilson, R.W., Dhillon, V.S., & HanifT, C.A. 1997, MNRAS, 
291, 819 



Bcrtrand Plcz and David L. Lambert: The Outer Atmosphere of Betelgeuse 



11 



Young, J.S., Baldwin, J.E., Boyscn, R.C., Haniff, C.A., 
Lawson, P.R., Mackay, CD., Pearson, D., Rogers, J., St.- 
Jacques, D., Warner, P.J., Wilson, D.M.A., & Wilson, 
R.W. 2000, MNRAS, 315, 635 

Young, K., Phillips, T.G., Knapp, G.R. 1993a, ApJS, 86, 517 

Young, K., Phillips, T.G., Knapp, G.R. 1993b, ApJ, 409, 725 

Weymann, R. 1962, ApJ, 136, 844 



